hep-ph/0509161 
TIFR/TH/05-35 

Event shape discrimination of supersymmetry from 
large extra dimensions at a linear collider 



Partha Konar* and Probir Roy 

Department of Theoretical Physics, Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumbai -400 005, India. 
E-mail: konar@theory.tifr. res. in, probir@theory.tifr. res. in 



Abstract 

The production of a charged lepton (i = e,fx) pair with a large missing energy at 
a linear collider is discussed as a means of distinguishing the minimal supersymmetry 
(MSSM) scenario from that with large extra dimensions (ADD) for parameter ranges 
where the total cross sections are comparable for both. Analyses in terms of event shape 
variables, specifically sphericity and thrust, are shown to enable a clear discrimination 
in this regard. 
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Introduction : A general expectation in high energy physics today is that of physics beyond 
the standard model (BSM) emerging at TeV energies. Supersymmetry (SUSY) [1] and 
extra dimensions [2] are two alternative possibilities in this direction that are the most 
exciting. They both address the naturalness/gauge hierarchy problem, arising from quantum 
corrections to the Higgs parameters, via the introduction of new physics at the TeV scale. 
Moreover, their attractive phenomenological features, in particular their promise of new 
states a bit beyond the current experimental lower mass bounds, put them in the limelight 
among scenarios of BSM physics to be explored by search strategies presently being designed. 
The latter, in fact, constitute the major motivation for constructing the next generation of 
colliders. If either SUSY or an extra dimensional scenario should manifest itself at sub- 
TeV to TeV energies, its signals ought to show up at the upcoming Large Hadron Collider 
(LHC) at CERN. It is widely accepted, nonetheless, that the precise nature of the BSM 
physics responsible for such signals may not always be easily gleaned from analyses of the 
corresponding data on account of the complexity of the hadronic environment in any LHC 
process. Indeed, in order to unambiguously identify the nature and detailed properties of 
any such new physics, a high energy and high-luminosity e + e~ machine [3] - such as the 
proposed International Linear Collider (LLC) or the Compact Linear Collider (CLIC) - will 
be very useful. 

We consider the signal comprising unlike-sign dielectrons/dimuons, produced in a linear 
collider together with a very high amount of missing energy, seeking to distinguish be- 
tween SUSY and and the Arkani-Hamed-Dimopoulos-Dvali (ADD) model [4] of large extra 
dimensions 1 ' 2 . Such a process has already been considered [7, 8] in the context of the Uni- 
versal Extra Dimension (UED) scenario [9]. The mechanism for this reaction is somewhat 
similar in SUSY and UED : a two-body production of heavy entities, each of which then has 
a dominant two-body decay. But the spins of the primarily produced entities are different in 
the two cases, leading to [8] distinguishable angular distributions and asymmetries. There 
are also differences in the lepton energy spectrum. We find, however that these quantities 
are not very sensitive to a SUSY vs ADD discrimination. First of all, the difference in these 
between the two scenarios is more quantitative, being in detailed shape aspects, rather than 
being something qualitative; systematic uncertainties would tend to wash out such quanti- 
tative differences. Secondly (and more importantly), these quantities are quite ISR-sensitive 
so that ISR-corrections significantly reduce the sensitivity to such a discrimination. 

1 Within the extra dimensional paradigm, there are other scenarios such as warped (Randall-Sundrum) 
or universal extra dimensions, which we do not address here. 

2 Another process where the two scenarios have been compared is e + e~ — > Jh$ where 7^ is a hard photon. 
The reactions for the ADD and SUSY scenarios in standard notation are e+e~ — ► G n -fH and e + e~ — ► GG^h 
[5] or e+e~ — ► XiXiIh [6] respectively. The energy spectrum of the hard photon together with the scaling 
of the cross section with CM energy and moment distributions of the transverse energy squared have been 
used for discrimination purposes. However, since there is only one observable particle in the final state, no 
event shape analysis is possible here. 
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Let us give an illustration to highlight the last point. The famous box-shaped lepton 
energy spectrum in the SUSY case has been found (as shown in Fig. 5 of [8]) to be squeezed 
in energy, looking more like a peak, after ISR corrections. When we compare this corrected 
spectrum with the peaked one for the ADD case, there does not seem a whole lot of difference. 
Similarly, the angular distributions are flat in either case for the bulk of the measurable range 
in the cosine of the angle between the two leptons. We do not include these plots here since 
that will detract from our central point which is the following. Distributions in event shape 
variables, such as sphericity and thrust, are known to be ISR-stable and are yet found to be 
sensitive to such a discrimination. They are qualitatively different between SUSY and ADD, 
having a peak in sphericity or break in thrust for the former and monotonic fall or rise for 
the latter. This is owing to differences in the mechanisms leading to the £ + £~$ final state 
in the two cases. Of course, slepton pair-production for SUSY will have a distinct threshold 
in ^/s unlike the generation of the corresponding ADD final state, the cross-section for which 
increases smoothly with y/s. So a careful scan of the CM energy for a threshold will also 
help discriminate between the two. However, that will require a more detailed step-by-step 
analysis. It will be useful to have a discriminant just with the first set of data at a particular 
\fs (above the slepton pair production threshold) and this is what we provide. 

We work within the minimal weak-scale R-parity conserving supersymmetric standard 
model (MSSM) which predicts the pair-production of charged sleptons [10], once the requisite 
energy threshold is reached, in an e + e~ collider. Each produced slepton would perforce decay 
into a charged lepton and the lightest supersymmetric particle (LSP). The latter is normally 
taken to be the lightest neutralino \i which, being stable and interacting only weakly, escapes 
unobserved through the detector - carrying a considerable amount of missing energy. In 
contrast, the ADD model has d extra dimensions compactified on a ci-torus. Together with 
time and the three spatial dimensions of our world, these constitute the bulk spacetime. 
The radius 3 R c of compactificaion of the extra dimensions could be as large as a quarter of 
a millimeter [11]. However, the SM fields are confined to a thin (thickness not more than 
1CT 17 cm [12]) .D 3 -brane, which is a soliton solution of the underlying string theory on which 
the ends of open strings are confined. A crucial feature of this model is that gravity, which is a 
property of spacetime itself, is free to propagate anywhere in the bulk. On compactification, 
a Kaluza-Klein tower of closely spaced gravitons appear in our spacetime, a large number 
of which (controlled by y/s) are producible 4 in a collider process [13] but are then invisibly 
lost in the higher dimensional bulk. To an observer on the brane, they would appear to be 
escaping unobserved with a large missing energy. This is a direct production of a three-body 
final state unlike the SUSY case where the decays of the heavy sleptons tend to generate 

3 For simplicity, we take the same radius of compactificaion for each of the d dimensions. 

4 An alternative way of probing the ADD scenario is to consider virtual graviton exchange [13] in SM 
processes where a coherent sum over closely spaced gravitons is involved, leading to deviations from SM 
predictions 
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more isotropic events. 

Comparison of the two signals: Recall that our process is e + e~ — > £ + £~f! where i sums 
over both e and /i. Charged slepton (e.L,R or Jil,r) pair production in an e + e~ collider with 
both unpolarised and polarised beams has been explored earlier [10]. Once produced, the 
sleptons decay into either a chargino-neutrino pair or into a neutralino-lepton pair. The 
partial decay widths are governed by both the mass and the composition of the charginos 
(neutralinos) as well as by the type (L or R) of slepton. We select the channels yielding the 
final state of a same-flavour unlike-sign dilepton associated with a missing energy 5 , namely 

e + e" - £ + L J- L , R - ttxlxl (1) 

In our analysis, we do not adhere to any particular SUSY-breaking scenario and make no 
assumption related to any high scale physics other than adopting gauge coupling unification. 
Thus, whereas the slepton masses 6 are free parameters in our analysis, the neutralino 
masses and couplings are completely specified by the respective SU(2) and U(l) gaugino 
masses M 2 and M 1 , the Higgsino mass parameter \i and tan/5, which is the ratio [1] of the 
two Higgs vacuum expectation values arising in the MSSM. 

The branching ratio for slepton decay into the lightest neutralino and the corresponding 
lepton depends on quite a few parameters: m^, //, tan (3 as well as the gaugino mass param- 
eters Mi and M 2 . Of these, the dependence on tan/5 is the least pronounced and therefore 
we shall henceforth use only one value of it, namely, 10. Thus, only four parameters remain, 
namely fi, M 2 and Mi. For a given slepton mass, the relevant branching fraction is then 
governed essentially by two factors: (i) the composition of the LSP and (ii) the energywise 
accessibility of slepton decay channels into the heavier neutralinos/charginos. The resulting 
dependence is still quite intricate and can be followed from ref [14]. 

Turning to the ADD scenario, the production of a dielectron or dimuon pair with missing 
energy has been considered [15, 16] earlier for probing its parameter space. The relevant 
reactions are 

e + e - ^ £+£-G n , (2) 

where n is summed incoherently over the energywise accessible part of the tower of closely 
spaced gravitonic modes. Two parameters, determining the relevant cross-sections, are: (i) 
the number of extra dimensions d and (ii) Planck's constant in the bulk or the so-called 
higher dimensional string scale Ms, expected to be in the TeV range. 

The SM backgrounds to our signal (of a same-flavour, unlike-sign dilepton pair plus a 
substantial amount of missing energy) arise from all processes of the form 

e+e" -> tt'ViVi (3) 

5 In case > M%±, there is also the chain e+e - — > ^lr^lr ~^ X v X~ v t^l ~^ ^^"XiXi^^z 7 ^- 
However, it makes a very small contribution, which we do take into account. 
6 Again, for simplicity, we take — m i — m i ■ 
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Table 1: Cross-sections (in fb), after the imposition of the cuts described in the text, for the 
e + e~ — > £ + £~$ signal in both ADD and SUSY scenarios at a linear collider with y/s = 500 
GeV (upper panel) and y/s = 3 TeV (lower panel). Parameter values are given in bold. 

where i can be any flavor. A significant portion of this background originates from the ££Z 
final states, with the real Z boson decaying into a neutrino pair carrying missing energy. 
These can be easily removed by imposing a suitable cut on the vv (missing) invariant mass. 
On the other hand, the background from e + e~ — > W + W~ £ + £~U£U£ can be explicitly 
subtracted by reconstructing (with a two-fold ambiguity) events with the on-shell W-pair. 
Such a procedure is, of course, problematic if there is an accompanying photon or if one of 
the W's is off-shell. But then the use of appropriate longitudinally polarised beams would 
lead to a drastic reduction of this type of background. Though we perform our present 
analysis with unpolarised beams, we shall comment on the use of beam polarisation at the 
end. 

Let us consider an e + e~ collider operating at a centre-of-mass energy 7 of 500 GeV (3 TeV). 
The kinematic cuts used in our analysis are as follows: 

• Each of the final state charged leptons should be at least 10° away from the beam pipe. 
This tames collinear singularities arising from t-channel photon exchange. At the same 
time, the elimination of any background effect from beamstrahlung is mostly ensured 

7 We here consider phase 1 of the ILC with CLIC taken up within brackets. 
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[17]. 

• Each charged lepton should have a transverse momentum p l T > 10 GeV (20 GeV). 

• We demand a missing transverse momentum p™ lss > 15 GeV (25 GeV) since this also 
helps 8 in reducing the two-photon background. 

• The tracks of the two unlike-sign leptons must be well-separated, with AR > 0.2, 
where AR = ^Arj 2 + A<p 2 in terms of the differences of the lepton pseudorapidities 
At/ and the azimuthal angles A<p. 

• The opening angle between the lepton tracks is required to be limited by 5° < 9 f + f - < 
175°. This ensures not only a sufficient missing energy, but also the elimination of 
possible cosmic ray backgrounds. 

• The missing invariant mass M miss has to satisfy the inequality |M miss — Mz\ > 10 GeV in 
order to eliminate the background of two final state neutrinos arising from Z-decay. A 
further increase in the value of the lower cut in M miss will reduce other SM backgrounds 
(such as those from off-shell Z and W + W + production). On the other hand, too high 
a value of M miss will reduce the signal. We have chosen a lower cut of 150 GeV (450 
GeV) on M miss which should reduce the SM background and at the same time yield a 
reasonable signal. 

The total cross-sections for our signal in the SUSY and ADD cases are presented in Table 1 
(with sample choices of parameters) after applying the event selection criteria described 
above. We see that these can be comparable in magnitude for sizable parametric regions. 
For the SUSY case, we have checked that our selectron pair-production cross sections match 
with what are plotted by Freitas et al. [10] in their Fig. 3(a) and moreover that the branching 
fractions for the decays of charged sleptons into charged leptons plus LSP, used as multiplying 
factors, agree with those of Choudhury et al. [14]. A similar check for the ADD case has 
been made with the results of [15]. The total contribution from SM backgrounds after 
applying the same event selection criteria comes out to be 36.4 fb (72 fb). One can estimate 
the minimum value of the signal strength for which significance, defined by S/ VB, where 
S (B) = number of signal (background) events, takes the desired value of 3. Considering an 
integrated luminosity of 100 fb -1 (1000 fb -1 ), the minimum signal strength needed for our 
analysis is found to be as — 1.8 fb (0.8 fb). 

Event shape variables: Though event shape analyses have so far been fruitful with many 
particle final states, even for two visible final state leptons, associated with a large missing 
energy, we expect to utilize the same to distinguish the signal of supersymmetry from the 
one originating from large extra dimensions of the ADD model. The underlying idea is 

8 cf. § III of rcf. [15] 
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Figure 1: Distributions of sphericity S (upper panels) and thrust T (lower panels) for the 
e + e~ — > C. + C.~]p signal at y 7 ? = 500 GeV both in the ADD (left panels) and SUSY (right 
panels) scenarios. The distributions are shown for different choices of parameters (M 2 , M 1; 
H, mi) in GeV for SUSY and (Ms, d) in the ADD case. 

the following. In the supersymmetric scenario a pair of rather heavy sleptons (£ + £~) are 
produced not far from the threshold in the center of mass energy of the collider. Each of 
these then subsequently decays into a lepton and a LSP. Because of the lack of any significant 
boost for each slepton, the daughter leptons lead to more isotropic events in the CM frame. 
The heavier the slepton, as compared to the LSP mass, the more isotropic is the event. In 
the ADD scenario, on the other hand, the two leptons are produced in association with a 
single graviton; in fact, a significant part of the cross-section comes from the production of 
heavier graviton modes which are energywise accessible. Hence, the events here are more 
spiked. It should be noted further that the event shape would also carry some signature of 
the spin information of the graviton. 

Event shape characteristics can be specified by constructing variables like sphericity, 
thrust, oblateness [18] and circularity [19]. It turns out to be sufficient for us to consider 
the first two. Criteria, made with such variables, have earlier proved to be promising, for 
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instance, in reducing multiple jet backgrounds to the tt signal at the Tevatron [19]. Sphericity 
is constructed from a normalised tensor, 

qij _ J2a PgPj ( A_\ 

where, p l a is the ith component of the three-momentum p a of the ath visible final state 
particle, with i,j = 1, 2, 3 and a summed sum over all such particles. S 1 ^ has three eigenvalues 
^1,2,3 with Ai > A 2 > A 3 > and Ai + A 2 + A 3 = 1. Since only two of the A's are independent, 
the sphericity of an event can be defined as S = §(A 2 + A 3 ). This is essentially a measure 
of the summed square of transverse momenta with respect to the event axis, two extreme 
cases being S — 1 for for an ideal spherical event and S = for a linear event. On the other 
hand, thrust is defined by the quantity, 

where \n\=l and this vector n is the thrust axis for which maximum is attained. The allowed 
range is | < T < 1, where a spiked shape event has T ~ 1 and an isotropic event corresponds 
toT~i 

Results and Discussions : We now analyse the sphericity (S) and thrust (T) distributions 
for our signal computed with tree level diagrams. The computation has been performed 
with the COMPHEP [20] program package in the supersymmetric case and with the HELAS 
subroutine [21] in the ADD case. In Figs. 1 and 2 we present these distributions for different 
choices of parameters. First, we note that, since there are only two visible particles here, the 
events are planer (A 3 = 0) with the shape being circular rather than spherical for 9 S max = |. 
The sphericity distributions in the case of SUSY are seen to clearly depend on the slepton 
mass rrig, events become more and more circular for larger and larger values of mi, showing 
a peaked structure, the peak location shifting towards the right as mi increases. On the 
other hand, sphericity distributions in the ADD case show a strong maximum at S — 0, 
monotonically falling faster with S, being largely insensitive to the values of Ms and d. 
We find in the SUSY case also that the sphericity distribution which is controlled 10 by the 
slepton mass m^, is insensitive to variations in M 2 ,Mi,/i. Coming to thrust distributions, 
we see that they peak at T = 1 for the ADD case and again do not change much if M s 
and d are varied. In the SUSY case, the thrust distributions are flatter except when mi is 
close to the present lower bound in which case they tend to resemble the ADD distributions. 
It is noteworthy that, even for as low as 110 GeV, the sphericity peak in the SUSY 
case is distinguishable from the maximum at S — for ADD. It is, however, true that the 
thrust distributions do not yield any additional advantage over the sphericity ones in so far 

9 With the ISR/FSR on, S can go all the way to 1. 

10 Thcrcforc, the value of the slepton mass rug can be extracted from a measurement of the location of the 
sphericity peak. 
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Figure 2: Distributions of sphericity S (upper panels) and thrust T (lower panels) for the 
e + e~ — > £ + £~$ signal at yfs = 3 TeV both in the ADD (left panels) and SUSY (right 
panels) scenarios. The distributions are shown for different choices of parameters (M 2 , M 1; 
H, mi) in GeV for SUSY and (Ms, d) in the ADD case. 

as the basic discrimination is concerned. We have also looked at oblateness and circularity 
distributions which show similar features with no additional advantage. We have checked 
that similar features characterize the \fs = 3 TeV case. 

One shortcoming of the above analysis may be that it has not included corrections due to 
the emission of collinear photons from the initial state (ISR) as well as the final state (FSR). 
Any sphericity distribution is regarded as vulnerable to changes caused by FSR because of 
the quadratic form of S^. Such an instability may be less pronounced for thrust which is 
linear in the concerned momenta. Though these effects are expected to be much smaller 
in leptonic processes than in hadronic ones because of the smallness of olem i n comparison 
with ctQCD, it is a concern that cannot be overlooked. We have, therefore, analysed these 
distributions for the SUSY case using the event generator pythia [22] which has the provision 
of including ISR and FSR contributions. The results, (with and without ISR and FSR) are 
shown in Fig. 3 for the y/s = 500 GeV case. Here the distribution of events, normalised by 
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Figure 3: Sphericity S (left panel) and thrust T (right panel) distributions for the e + e~ — > 
£ + £~ ]fi signal at ^ = 500 GeV in the SUSY case reproduced in pythia to include ISR and 
FSR effects (thin lines are without ISR and FSR effects). The different choices of parameters 
in the distributions are as in Fig. 1 



the total number of events times the bin size 11 , are considered for sphericity and thrust. We 
notice some small but visible changes in the sphericity distributions with and without ISR 
plus FSR - but not enough to adversely affect the discrimination from the ADD case. As 
expected, all such changes are significantly less in the thrust distributions. 

A most important aspect of Fig. 3 is the following fact: unlike the shape of the lepton 
energy spectrum [8], the locations of the sphericity peaks for the SUSY case are unaffected 
by ISR and FSR effects. Moreover these locations are the same 12 in Figs. 1 and 3. This 
location is therefore a robust discriminant between the SUSY and ADD cases since, for the 
latter, the peak is always at S = 0. Fig. 4 shows a correlation plot for ^Js = 500 GeV 
between the total cross-sections for the process e + e~ — > £ + £~ $ and the locations of the 
sphericity peaks shown as scatter points for some sample parametric choices. While the cross 
sections for the SUSY and ADD scenarios overlap, the sphericity peaks are very distinctly 
apart. The maximum for the ADD case is at S = 0, while for SUSY it is at a nonzero value 
of 5* which is an increasing function of m^- moving towards the limit | when approaches 
^Y- We have rescaled this whole analysis by taking yfs = 1 TeV with higher values of the 
concerned parameters and have found very similar results. Thus even if SUSY or a model 
of extra dimensions a la ADD is beyond the reach of the first phase of the ILC and is only 
accessible to its second phase or to CLIC, our method of discrimination should work. 

11 We have chosen the same bin size for Fig. 1 and Fig. 3 

12 The agreement between the plots of Fig. 1 and the thin lines of Fig. 3 is a cross-check on the consistency 
between the calculations with COMPHEP and with pythia. 
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As mentioned earlier, the use of appropriate longitudinally polarised beams would sup- 
press the WW background significantly and one could have a reach for larger parameter 
spaces in both the SUSY and ADD scenarios. However, the distinctive features of sphericity 
and thrust distributions, as noted here, would be unaffected since their origin has nothing 
to do with beam polarisation. Let us also comment on beamsstrahlung [23] effects which 
we have ignored. Because of our lower cut on p™ ss , beamsstrahlung photons going down 
the beamlines will not affect our analysis. Nevertheless, beamsstrahlung will cause a degra- 
dation in the effective value of the CM energy ^Js. This will quantitatively affect precision 
measurements, such as the extraction of mi from the location of the sphericity peak, causing 
a systenmatic error which will need to be taken into account. However, from the estimates 
made in [23], we deem it unlikely that this degradation will change the qualitative differ- 
ence in the sphericity distribution between the SUSY and ADD scenarios, namely a peak 
structure in the former and a strutureless monotonic falloff in the latter. 
Summary and Conclusions: In the location of the maximum in the sphericity distribution 
of the process e + e~ — > £ + £~f!, we have discovered a sharp and robust discriminant between 
the SUSY with ADD scenarios which should be utilisable either at ILC or at CLIC wherever 
these scenarios become accessible. Longitudinally polarised beams with appropriate helici- 
ties should help further by suppressing the SM background. Similar considerations can be 
extended to the subprocess [24] qq — > £ + £~pr , accessible at the LHC. However, a careful 
calculation of QCD corrections to it, with particular reference to ISR gluons, will be required 
first. 
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